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Abstract. We provide details for the proof of convergence of the locally inertial 
Godunov method with dynamical time dilation^ introduced by Vogler in |13) . 



1. Introduction 

In this paper we fill in the details in the proof of convergence stated in Section 7 of 
for the locally inertial Godunov method with dynamic time dilation, a numerical 
method for computing shock wave solutions of the Einstein equations in Standard 
Schwarzschild Coordinates (SSC). We refer the reader to [II] for an introduction and 
for the notation assumed at the start here. 

The main conclusion of the theorem stated and proved below is that a sequence of 
approximate solutions (max, Aax) — ^ (m, A) of the locally inertial Godunov scheme 
that converge boundedly without oscillation to a limit function (m,A), must be a 
weak solution of the Einstein Euler equations in SSC, c.f. [HI [2], 

ut + /(A, u)^ = g{A, u,x), 

A' = h{h,u,x). ^ ■ ' 

That is, the theorem reduced the proof that one has an exact solution to the two things 
most easily established in a numerical simulation: convergence without oscillations. 

System (11.11) is the locally inertial formulation introduced in [2], which is weakly 
equivalent to the Einstein equations, (see [Hj for discussion). The proof here is a 
modification of the Groah and Temple argument used for the locally inertial Glimm 
scheme [2] , with several differences. Both the Glimm scheme and the Godunov meth- 
ods employ Riemann problem approximations, but the main difference is that the 
Godunov method employs averaging rather than random sampling at the end of each 
time step. The theorem proved here assumes a total variation bound, (it is still an 
open problem to prove such a bound), while the Groah and Temple theorem estab- 
lishes this bound by an argument using wave strengths to bound the total variation of 
waves in the Riemann problem solutions. In the theorem here, we also allow variable 
time steps. The final difference is that here we include right and left boundary data 
reflecting the limited extent in space of a computer simulation. 



This work summarizes resuhs credited to Vogler's doctoral dissertation, (13) . which was supervised 
by Blake Temple. Both authors were partially supported by second author's NSF Grant, where the 
problem was first proposed. 
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There are two main steps in the proof. The first is to show the discontinuities in 
the metric A along the boundary of Riemann cells are accounted for by the inclusion 
of the term 

A'- Va/(A,,,m,x) (1.2) 

in the ODE step ^14j. The second step is to prove the jump in the approximate 
solution along the time steps are of order Ax. In their work Groah and 
Temple did not need the convergence and total variation assumptions because with 
the Glimm scheme, these assumptions are proven as long as there exists a total 
variation bound on the initial data, a truly remarkable feature of the scheme. In 
this paper, applicable to the simulations in fl3\ [H], these assumptions are natural 
assumptions that can be verified numerically. In particular, the theorem is perfectly 
suited to the numerical simulation of points of shock wave interaction in fiSi [H]. 
Once one numerically establish convergence and a total variation bound, the theorem 
here implies convergence to a weak solution of the Einstein equations. 

2. The Convergence Theorem 
The main theorem of this paper is the following: 

Theorem 2.1. Let UAx{t,x) and Aj\x{t,x) be the approximate solution generated 
by the locally inertial Godunov method starting from the initial data UAxito,x) and 
AAx{to,x) for to > 0. Assume these approximate solutions exist up to some time 
tend > to and converge to a solution {uax, Aax) — ^ {u, A) as Ax — > along with a 
total variation bound at each time step tj 

T.V.lr^^„,r^.A^Ax{t„ ■)}<¥, (2.1) 

where T.V.i^r^^„,rmax]WAx(tj, ■)} represents the total variation of the function UAx(tj,x) 
on the interval [r min i"^ max]- Assume the total variation is independent of the time step 
tj and the mesh length Ax. Then the solution {u, A) is a weak solution to the Einstein 
equations (1.26)-(1.29) in |2]. 

Proof. Suppose we have approximate solutions (max, Aax) obtained by the locally 
inertial Godunov method that satisfy the hypothesis of the theorem. Having a total 
variation bound at each time tj places a total variation bound on the inputs to all 
the Riemann problems posed at that time. In [2], Groah and Temple show a total 
variation bound on the inputs implies a total variation bound on the solution to the 
Riemann problem for any time t such that tj < t < tj+i. By the self similarity of 
the solution to the Riemann problem, this result also implies a total variation bound 
for any space coordinate within the Riemann cell. More specifically, we have the 
following bounds: 

r.K[.,.„.a«A.(t,-)}< v^, (2.2) 

and 

T.V.[t„t,^,){uAx{;x)}<V, (2.3) 

for any x and t within the Riemann cell Rij. 

All the functions /, G, and g derived in [2] are smooth, and it is the metric that 
is only Lipschitz continuous. The smoothness of these functions is used throughout 
this proof. 
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Let T = tgnd — to be the overall time of the solution, and for each mesh length Ax 
define the minimum time length 



At = mm{Atj} (2.4) 

j 



as the minimum over all the time lengths defined in p3:j. By definition, this time 
length is proportional to the mesh length. At oc Ax, implying 0{At) = 0{Ax), 
and there exists a constant C bounding all the time lengths, Atj < CAt for all j. 
Throughout this section, let C be a generic constant only depending on the bounds 
for the solution [to,^end] x [rmin^fmax]- This variable is created to unify all the time 
steps, and more importantly, used to calculate the maximum number of time steps 
needed to go from to to tend- 

We now follow the development of Groah and Temple in [2]. Recall from [Ti] . 
M^^(t,x) denotes the collection of the exact solutions in all the Riemann cells Rij for 
the Riemann problem of the homogenous system 

ut + f{A,,,u)^ = Q. (2.5) 

So M^^(t, x) satisfies the weak form of this conservation law in each Riemann cell 

^= {-u^x^t - /(Ay, u^^)ipx} dxdt 

J J Rij 

JR^ (2.6) 

+ / {f\Aij,ul^xit^Xi))vit^Xi) 
jRj 

-/(Aij, (t, Xi-i))(p{t, Xi-i) } dt, 
where is a smooth test function with Supp{ip) C [to, tend) x [o, &] for a < Tmin < 

Remember from [T3], u{t,UQ) denotes the solution to the ODE 

ut = G{Aij, u, x) = g{Aij, it, x) - A' ■ VA/(Ay, u, x), 

U{0) = Uq. 

Therefore, 

u{t,uo) = Uo+ / {g{Aij,u{^,uo),x) - A' ■VAf{A,j,u{^,uo),x)}dC (2.8) 

^0 

Also, recall from [H] u/\x denotes the approximate solution obtained using the frac- 
tional step method. Since our fractional method takes the Riemann problem solution 
and feeds it into the ODE step, max is defined on every Riemann cell Rij as 

UAx{t,x) = uf^{t,x) + / {g{Aij,u{^ -tj,u^^{t,x)),x) 

(2.9) 

df 

This expression implies the error between the approximate solution and the Riemann 
problem solution is on the order of Aa;; a fact that is repeatedly used throughout the 
proof. 
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Define the residual e = e{uj\x, ^Ax, v) of "^Az and Aa^ as the error of the solution 
in satisfying the weak form of the conservation law f 1 1.1 1) by 

£(mAx, Aax, V^) = / / {-UAx'-Pt- f{A^^,UAx)'-Px - g{AAx,UAx,x)(p}dxdt 

-h-h 

i=n+l » 

/ {-UAx^t - fi-^ij, uax)^x - g{-^ij, uax, x)ip} dxdt 

: 1 ^ Rij 

(2.10) 



where 



n+l 



h= UAxito , x)dx = ^ UAxito^x)dx, (2.11) 

and 

h= {f{Aij,UAxit,r:^iJ)ifi{t,r^iJ - /(A^^, MAx(t, r+„^))<^(t, r+<,^)} dt 

J to 

= Y1 / {f(^ij^'^^^(t,r^J)ip{t,r^iJ - f{Aij,UAx{t,r^^J)ip{t,r^^j}dt, 

(2.12) 



The expression X]i=it^ denotes a double sum where the index i runs across all the 



spatial gridpoints, and the index j runs across all the temporal gridpoints. Recall 
from [Hj, n is the number of spatial gridpoints, and there are n + l Riemann cells. 
Our goal is to show e{uAx, Aax, V^) = 0{Ax) because if the approximation converges 
{uax, Aax) {u, A) as Ax — )■ 0, then the limit function satisfies the condition of 
being a weak solution to the Einstein equations e{u, A, (p) = 0. 
Substituting (^M> into d^TD]) gives us 

i=n+l „ „ 

£ = / / {-U^xft - f{Aij,UAx)fx - g{Aij,UAx,x)(p 



Rij 



-Vt [ [9{Aij,u{^-tj,u^'^{t,x)),x) ^ ■ ^ 

df 

- ^ (^^i' - h^^Axit, x))) ■ a'aJ d^} dxdt -h- h. 



Define 



Jtj 

df 



(2.14) 
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Plugging the weak form of the conservation law (12 .Gp of each grid rectangle into fl2.13p 
gives us 



1=71 + 1 



— iL-r-L p p 

^ {'^x[f{A.ij,uf^)-f{Aij,UAx)]-g{-^ij,UAx,x)(p 

i=i,j -^R^j 

-(ptlljit,x)} dxdt 

i=n+l „ 

i=l,j -^R^ 
i=n+l „ 

-h-^ / {f{Aij,u^^^{t,Xi))^{t,Xi)-f{Aij,u^l{t,Xi^i))^{t,Xi_i)]dt. 

(2.15) 



Note 



fiAij,uZ) - f{A,„UAx) \ < CAt 



(2.16) 



which implies 

i=n+l 



ijy'^Ax) ~ f{Aij^UAx)]dxdt 



<CM^AeAx(^\ (n + l) = 0(Ax 



(2.17) 



where the number of time steps is proportional to T/At and the number of space 
steps is 0(1/ Ax) pj. 

Since (t^,x) = UAx{t'j',x), the following sum is rearranged to become 



i=n+l 



-A- E / {'"Ar(^7+i.a;)v?(tj+i,a;) -u^^(t+,x)(^(tj,a;)}c?x 
{uAxitj,x) - u^^{tj,x)} Lp{tj,x)dx 

n 

<f{tj,x) {uAx{tj,x) - UAxit~,x)]dx 

n 

■,ax 

(p{tj,x) {uAxitJ^x) -UAx{tJ,x)}dx, 



(2.18) 



where the term UAxitj, x) is added and subtracted to isolate the jump in the solution 
uax across the time step tj. We define this jump ei = ei{uAx, Aax, V^) as 



£i{uax,Aax,v) = E 



(p{tj,x) {uAx(t^,x) - UAx{tj ,x)] dx, (2.19) 
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and this definition allows us to rewrite fl2.15p as 



i=n+l 



e=0{/\x)+ei+ ^ {-g{Aij,UAx,x)(p - (ptlljit,x)} dxdt 



j^O min 

i=n~\'l 



/ {/(Aii,MAx(^,a;i))(^(t,x,) - /(Aij,M2^(t,Xi_i))<^(t,Xi_i)}c/t 

(2.20) 



But the last sum is rearranged to cancel the boundary conditions as follows: 



i=n+l 



"^2- / {/(Aij,MAr(^,a;i))(/?(t,Xi) - /(Aij,M2^(t,Xi_i))(/?(t,Xi_i)}rft 

= 5Z / {/(Ai+i,i,«Ar(^.a;i)) - /(Aij,u2^(t,Xi))}(^(t,a;i)c?t 

+ / {/(Ai,i'«Ar(^>a;o)) - /(Aij,MAz(i,a;o))}</'(t,a;o)c/t 

+ 5]! / {/(^n+ij'^Ari^^a^n+i)) - /(A„+ij,MAx(t,a;n+i))} V5(t,a;„+i)c/t, 



(2.21) 



where 



/ {/(Ai,j,MAr(^)a;o)) - /(Aij,MAx(i,a;o))} V5(t,Xo)(it 



<ll^llooCAt2 =0(Ax 



(2.22) 



and similarly 



/ {/(A„+ij,u2^(t,x„+i)) - /(A„+ij,UAx(i,a;„+i))} V5(t,a:„+i)c?t = 0(Aa;). 

. -^^^ 

(2.23) 

Note that the resulting double sum in fl2.2ip lost a term, resulting in only n terms. 

To simplify the l\- term, we add and subtract a term deviating from it by an order 
of Ax, use integration by parts on the new term, and with the result add and subtract 
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another term to reduce the expression further. To this end, let 



i=n+l 

I AS = 



=n+l „ „ „t 
j=l,7 J Rij Jtj 

df 

df 

(^^i' ""(^ ~ ^' ""A^l^, x))) ■ ^Ax} didxdt. 



OA 

OA 



(2.24) 



From the total variation bound on the Riemann problems and the smoothness of /, 
this term is bounded by 



i=n+l 



\Ias\< ^ \Moo C T.V.[^._^^^.]{uAx{-,tj)}d^dxdt 

i=l,j Rij 

< CAt'Ax J2 tj)} (2.25) 



<CVyt\LAxAt'^^=0(Ax'), 



and the above procedure reduces the term to 



i=n+l 



p p ^ — /fc-T-L p p pt 

- (ptllj{t,x)dxdt^ Ias - / / / {ai^ij^H^-'tj^f^Axi^^^))^^) 

J JRij i=l,j *J 



df 

- (A.„ M(e - tj,ul^,{C, x))) ■ A'^ J dxdt 



0{Ax')-Yl / [g{Aij,u{C-tj,ul^,{C,x)),x) 
df 

- - ^^^(e, a;))) • A'^,] di 

- / liAx, a;) - ^(Ay ' ""Ax) • A'^Jd^ 

i=n+l ( tj+i 

0(Aa;2) - J2 \ '^(^i+i' ^) / [^(A^^' - "A^l^.+i, a;)) 



9A 



Aij, «(C - tj, M^^(ti+i, x)))- A'^ d^^dt + h + h, 

(2.26) 
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where 



i=i,j -^R^ I -^h 



X) 



g{Aij,u{^ -tj,ul^{^,x)),x) - -g^iAij,u{^-tj,u^^{tj+i,x))) ■ A'^^ 



(2.27) 



and 



i=n+l 




df 

g{Aij,UAx,x) - -g^{Aij,UAx) ■ A'^^ 



dxdt. 



(2.28) 



Again by smoothness and the total variation bound, we have 



i=n+l 



I/4I < II^IL XI ^ T.1/.[^^_^,^^] {MAx(-,ti)} AxAt 



T 



< ||^||^CAxAt5^T.\/.[,_,,_] {nA.(-,t,)} = M\^CV/\x/\t— = 0(Ax) 



(2.29) 



Substituting fl2.2ip and f l2.26p into fl2.20l) along with using 02.91) as an identity leaves 
us with 



i=n+l 



e=0{Ax)+ei- ^11 (p^{Aij,UAx) ■ A'^^dxdt 



i=n „ 

+ XI / <^(^,a;i) {/(Ai+ij,u^^(t,Xi)) - /(Aij,uf^(t,Xi))}rft 



(2.30) 



The second sum represents the jump in the flux function /, resulting from the dis- 
continuities in the metric A, and the first sum is the addition to the ODE step ( 12. 7p 
specifically designed to cancel these jumps in the flux. 

To see how the cancelation works, we perform a Taylor expansion on the test 
function, and we add and subtract terms deviating by order Ax. The first sum in 
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fl2.30p is expanded as 



i=n+l 



OA 

i=n+l 

i=n+l f. !• f Ci e Qj 

'-Ax ~ ' dxdt 



^ I [ v{xi,t)^{Aij,UAx) ■ A'^^dxdt + 0{Ax) 



/ v{xi,t) I —(Ay, Max) -A'^ 
+ E / ^{x^,t) f l^iA,^,uZ)-Kx-^iA,,,uZixi,t))-A'Adxdt 

i=i,j JRj ^ 

Aax{x + tj), u^^{xi, t)) ■ A^^l dxdt 



OA 

i=n+l 



\^ f f^' df Ax 

+ 22 V{x^,t) / -!-{Aax{x + —,tj),u^^{xi,t)) ■ A'^Jxdt + 0{Ax) 

i=l,j JRj ^ 



(2.31) 



From the smoothness of /, each of the first three sums in equation fl2.3ip are 0(Ax) for 
the following reasons: the first sum is order Ax from the ODE step in the definition 
of the approximate solution uax (12. 9p . the second sum is order Ax^ by the total 
variation bound on solutions to the Riemann problems, and the third sum is order 
Ax by the Lipschitz continuity of the metric A. After these bounds are established, 
(I2.3ip reduces to 



j=i,j J J Ri] 

j=n+l 



~ /" f^^ df Ax 

J2 / Vi^i^t) -!-{Aax{x + —,tj),ul^{xi,t)) ■ A'^Jxdt + 0{Ax) 

i=i,j JRj J^^-i ^ 
i=n+i ^ r' df Ax 

22 / —{Aax{x + —,tj),u1l{xi,t))dxdt + 0{Ax 

i=ij Rj Jxi-i 

i=n+l 

I /. \r^/» RP/. \\ ^/» RP / . \\1t. .^/a 

iX 



J2 [ v{t,x,){f{A,+,^^,ul^^{t,x,+))-f{A,„uZ{t,x,+))}dt + 0{Aa 
i=i,j -^Rj 

(2.32) 
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Plugging this result f l2.32p into fl2.30p gives us 
e = 0{Ax) +ei 

^ii,Xn+l) {/(A„+2,i,MAr(^>a;„+i)) - f{An+ij,U^^{t,Xn+l))} dt, 

(2.33) 

where one term remains due to the mismatch in the number of terms in the spatial 
sum. Clearly, this last term is 0(Aa;). 
So the residual boils down to 

e(MAx, Aax, ^) = £i{uax, Aax, ^) + 0{Ax), (2.34) 
with all that remains to show is 

£i = ^ / ip{tj,x) {uAx{t^ - UA^{t~ ,x)] dx = 0{Ax). (2.35) 
To estimate ei, we break up the sum by each time step tj and define 



fmax 



ip{tj,x) {uAx{t^,x) - UAx{tj ,x)}dx 

n 

/ V{t^,x){uAx{t],x)-UAx{t~,x)]dx, 



(2.36) 



with Xj+ = X,- , 1 and Xi- = . 

'"•"2 2 

Recall from [H] , the approximate solution for the new time step is computed by 
the Godunov step, using averages at the top of each Riemann cell Rij. In particular, 
the solution at each new time step is 

UAx{t^, x) = u{tj - tj^i,u{tj), x)) (2.37) 

where 

^(^^) = ^ f ^^uZ{tj,x)dx (2.38) 

To finish the proof, a lemma is needed, which is proven at the end of this section. 
This lemma states the difference of the ODE step taken on an average verses the 
solution to the Riemann problem across the top of the Riemann cell is bounded by 
the total variation of the Riemann problem. 

Lemma 2.1. Let represent the solution of the Riemann problem in the Riemann 
cell Ri,j-i and UAx(t) denote the average of the Riemann problem solution across 
Riemann cell. Let u be the solution obtained by the ODE step 1^2. 7| j and (f be a 
smooth test function. Then the following bound holds 



(2.39) 



{u{tj - tj_i,UAx{tj),x) - u{tj - tj_i, u^^^itj, x), x)]ip{tj, x)dx 

< C'll'^lloo^^^^ T.V.[x^,x,+^]{uAx{tj,-)} 

for some constant C. 
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Using Lemma I2.H (I2.36p is rewritten as solutions to the ODE step (12. 7p and 
bounded by 

(2.40) 



By the total variation bound on UAxitj, ■), the residual is bounded by 

At 



Bi < J^CIIv^lLAxAt T.K[,_,,_]K^(-,t,)} < C^^AxAtV = 0{Ax). (2.41) 



Therefore, e = 0{Ax) and the proof is complete. □ 

To prove Lemma 12. a preliminary result is needed: given a function on a set of 
points the difference of the function between any point and the average is bounded by 
the total variation of that function on the set. This result is provided by the following 

Lemma 2.2. Let u{x) be a function on the set [xi-,Xi^] and 

u = — u{x)dx (2.42) 



Ax 

he the average of u on this set. Then we have 

\u-u{x)\< sup \u{xi) -u{x2)\<T.V.^^^_^^.^^^{u{-)}. (2.43) 

Proof. The second inequality is true by the definition of the total variation 

sup \u{xr) - u{x,)\ < T.V.^,_,^^^{u{-)}. (2.44) 

X\,X2€[Xi,XiJ^l] 

To prove the first inequality, we assume it is false to obtain a contradiction, so suppose 
there exists x^ G such that 

sup \u{xi) — u{x2)\ < \u — u{x^)\. (2-45) 

X\,X2&[Xi,Xi + l\ 

Relabel the u-coordinates by an isometry ip : u v that maps the point u{x^) to 
the origin in the v-coordinates (i.e. (p{u{x^)) = 0), and the vector u — u{x^) in the 
direction of the 1st coordinate v^, as show in Figure [H 

Since the average of a collection of points is independent of the coordinate system 
in which they are labeled in, we have 

V = — — / v{x)dx = — — / (p{u{x))dx = (p I — — / u{x)dx j = (p{u) (2.46) 

J Xi- Jxi- \lS.X J / 

The following inequality holds by transforming equation f l2.45p over to v-coordinates 

\v{x)\ = \u{x) — u{x^)\ < \u — u{x^)\ = \v\ Vx G [xj_,Xj+], (2.47) 



which implies 
1 
Ax 



Xi + l 

v{x)dx 

Xl 



< / \v(x)\dx < — / \v\dx=\v\. (2.48) 
- Ax L_. ' ^ ^' Ax .L ' ' ' ' ^ ' 
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V2 



d 



Ui 



O = ^{u{x,)) 
Figure 1 . The isometry ^ : u v 



Vl 



This inequahty |f | < \v\ is an obvious contradiction, proving the first inequahty in 



dm. 



□ 



Now we prove the lemma used in the proof of Theorem 12.11 
Proof of Lemma 12.11 Recall the solution to the ODE step has the form: 

u{tj - u^^{tj, x),x) = u^^{t, x)+ G{Aij, u1l{t, x),x)dt. 



(2.49) 



This solution implies the LHS of (12.391) is written out as 

{u{tj -tj_i,U/^^(tj),x) -u{tj -tj_i,uf^{tj,x),x)}Lp(tj,x)dx 

+ I {G{Aij, UAx{t), x) - G{Aij, u^^{t, x), x))dt} ip{tj, x)dx 



(2.50) 



{uAxitj) - u^xi'tj^x)] 'f{tj,Xi)dx 
+ 1 / {G{Aij,UAxitj),x) - G{Aij,u^^{tj,x),x)} dt (p{tj,Xi)dx 

I Xi- -Jtj-l 

+ 0{Ax^), 

where the test function in the first term is approximated by a Taylor expansion. By 
the definition of the average function u, the first term is zero. By the smoothness of 
G, the bound (I2.39P is proven by 

{u{tj -tj^i,UAx{tj),x) -u{tj -tj^i,u2^{tj,x),x)}ip{tj,x)dx 



< C\\ip\\^AxAt sup {\uAxitj) - u^^{tj,x)\} 

< CyW^AxAt T.V.[^,_^^^^]{uAx{tj,-)}, 



(2.51) 
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where Lemma [2^2] is used to bound the difference between the average and the solution 
to the Riemann problem. 

□ 
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